Hyperfine-mediated transitions between a Zeeman split doublet in GaAs quantum 

dots: The role of the internal field 
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We consider the hyperfine-mediated transition rate between Zeeman split spin states of the lowest 
orbital level in a GaAs quantum dot. We separate the hyperfine Hamiltonian into a part which is 
diagonal in the orbital states and another one which mixes different orbitals. The diagonal part 
gives rise to an effective (internal) magnetic field which, in addition to an external magnetic field, 
determines the Zeeman splitting. Spin-flip transitions in the dots are induced by the orbital mixing 
part accompanied by an emission of a phonon. We evaluate the rate for different regimes of applied 
magnetic field and temperature. The rates we find are bigger that the spin-orbit related rates 
provided the external magnetic field is sufficiently low. 

PACS numbers: 71.23.La, 71.70.Jp 



I. INTRODUCTION 

The electron spin states in bulk semiconductor and 
heterostructures have attracted much attention in recent 
years. Experiments indicate very long spin decoherence 
times and-SHiall transition rates between states of dif- 
ferent spintxEltl These promising results have motivated 
proposals for information processing based on electron 
spins in quantum dots, Wjifich might lead to a realization 
of a quantum computerm! 

A quantum dot is region where electrons are con- 
fined. The energy spectrum is discrete, dueJ;g the small 
size, and can display atomic-like propertied!!'!! Here we 
will consider quantum dots in GaAs-AlGaAs heterostruc- 
tures. The main reasons for studying them are that rele- 
vant quantum dots are fabricated in such structures and 
Ga and As nuclei have a substantial hyperfine interac- 
tions with the conduction electrons. There are two main 
types of gate controllecLdots in these systems, so-called 
vertical and lateral dotso. They differ in transverse con- 
finement, which is approximately a triangular well and a 
square well for the lateral and vertical dots, respectively. 

Manipulation the electron spin while maintaining 
phase coherence requires that it should be relatively well 
isolated from the environment. Coupling a quantum dot, 
or any closed quantum system, to its environment can 
cause decoherence and dissipation. One of the measures 
of the strength of the coupling to the environment are the 
transition rates, or inverse lifetimes, between the quan- 
tum dot states. In GaAs there are two main mechanisms 
that can cause finite lifetimes of spin states. These are 
the spin-orbit interaction and the hyperfine interaction 
with the surrounding nuclei. If a magnetic field is applied 
the change in Zeeman energy accompanying a spin-flip 
has to be absorbed by phonon emission. For two elec- 
tron quantum dots, where the transitions are between 
triplet and singlet spin states, both spin-orbitfru and 
hyperfineEy mediated transitions have been studied. In 
both cases the transition rates are much smaller than the 
usual phonon rates, the spin-orbit rate being the highjsc 
except when the excited singlet and triplet states cross. til 



For low magnetic fields, i.e. away from the single-triplet 
transition, the energy of the emitted phonon can be quite 
large and the transitions involve deformation phonons 
rather than piezoelectric ones. 

If there is odd number of electrons in the dot the 
ground-state is usually a spin doublet so that energy 
change associated with the spin-flip is the electron Zee- 
man energy. Owing to the small g-factor in GaAs this 
energy is rather small compared to the level spacing and 
the dominating phonon mechanism is due to piezoelectric 
phonons. Recently spin-orbit mediated spin-flip tran- 
sitions between Zeeman levels was investigated!!!! Due 
to Kramer's degeneracy the transition amplitude for the 
spin-flip is proportional to the Zeeman splitting. This 
results in a spin-orbit spin-flip rate proportional to the 
fifth power of the Zeeman splitting. 

In this publication we consider hyperfine mediated 
transitions between Zeeman split levels in a quantum 
dot. The transition amplitude remains finite at zero ex- 
ternal magnetic field, giving a spin-flip rate (Eq. j27| ) 
) that is proportional to the cube of the Zeeman split- 
ting. The cause of this is an internal magnetic field due 
to the hyperfine interaction. We consider the impor- 
tant concept of internal magnetic field in some detail. 
Since the parameters of hyperfine interaction between 
conduction band electrons and undetbzicg nuclei in GaAs 
have been extensively investigatecuJ'Lj, including the 
Overhauser effe<^p«i,pepr relaxation in GaAs/AlGaAs 
heterostructures^taE3l!3l!3 , we are able now to predict 
the typical time scale for this process in particular quan- 
tum dot geometries. 

Upon completion of this work wc learned about recent 
results of Khaetskii, Loss and GlazmannO. They con- 
sider essentially the same situation and model and obtain 
electron spin decoherence without considering any mech- 
anism of dissipation. This is in clear distinction from the 
present result for spin-flip rate that requires a mechanism 
of dissipation, i.e. phonons. 

The rest of the paper is organized as follows: in sec- 
tion II the model used is introduced in addition to the 
basic assumptions and approximations used, section III 
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deals with the internal ma gnet ic field due to the hyper- 
fine interaction and section [V contains the derivation of 
the transition rates. Finally, in section [v| the results are 
discussed. 



II. MODEL AND ASSUMPTIONS 

We consider a quantum dot embedded in a Al- 
GaAs/GaAs heterostructure. The details of the quantum 
dot eigenstates are not important for now, it suffices to 
say that the energy spectrum is discrete and the wave 
functions are localized in space. The spatial extension 
of the wavefunction in the lateral and transverse direc- 
tion (growth direction) are denoted with i and z Q , re- 
spectively. Quantum dots in these heterostructures are 
formed at a GaAs/AlGaAs interface, where the confining 
potential is very strong so that t ;§> zq. We define the 
volume occupied by an electron as Vqd = tt£ 2 zq. The 
Hamiltonian of the quantum dot can be written on the 
form 



H = J2(ei+9»BB-S)\l){l\, 



(1) 



where £i are the eigenenergies which depend on the struc- 
ture of the confining potential and the applied magnetic 
field. The magnetic field also couples to the electron spin 
via the Zeeman term. 

Since the T point of the conduction band in GaAs is 
mainly composed of s orbitals the dipole interaction with 
the nuclei vanishes and the hyperfine interaction can be 
described by the usual contact term 



H HF = AS-J2lkS{r-R k ), 



(2) 



where S (Ik) and r (Rk) denotes, respectively, the spin 
and position of the electron (fcth nuclei). The delta func- 
tion indicates the point-like nature of the contact inter- 
action will result in a position dependant coupling. The 
coupling constant A has the dimension Volume x Energy. 
To get a notion of the related energy scale, it is straight- 
forward to relate A to the energy splitting of the doublet 
for a fully polarized nuclear system, 



E n — AC n I, 



(3) 



C n being density of nuclei and I is the spin of a nucleus. 
In GaAs this energy is E n w 0.135 me Vp-, which corre- 
sponds to a magnetic field of about 5T£j For a given 
quantum dot geometry the number of nuclei occupying 
the dot is defined as -/Vqd = C„Vqd- Since the Ga and 
As nuclei have the same spin and their coupling constants 
are comparable, we will assume that all the nuclear sites 
are characterized by the same hyperfine coupling A, and 
the C n w <Iq 3 where a is the lattice spacing. For real- 
istic quantum dots -/Vqd ~ 10 4 — 10 6 , and it is therefore 
an important big parameter in the problem. 



The coupling between the electron and the phonon 
bath is represented with 



H l ^ = Y j aM){K^ lqr 



(4) 



where w qv and b qL , are creation and annihilation operator 
for the phonon mode with wave-vector q on brach v. In 
GaAs there are two different coupling mechanisms, de- 
formation ones and piezoelectric ones. For transitions 
between Zeeman split levels in GaAs, i.e. low energy 
emission, the most effective phonon mechanism is due 
to piezoelectric phonons. We will assume that the het- 
erostructure is grown in the [100] direction. This is the 
case for almost all dots and it imposes important sym- 
metry relations on the coupling coefficient. The square 
of the coupling coefficient for the piezoelectric phonons 
is then given by (see Ref. |2l|) 



al{q) 



(ehi^h 
2pc w Vq 



AM 



(5) 



where (e/114) is the piezoelectric coefficient, p is the 
mass density, c v is the speed of sound of branch v, 
V is normalization volume and we have defined q = 
g(cos (f>s'm 9, sin c/>sm9, cos 9) and the A v 's are the so- 
called anisotropy functions, see appendix 0. 



III. INTERNAL MAGNETIC FIELD 

In this section we will introduce the concept of the ef- 
fective magnetic field, the internal field, acting on the 
electron due to the hyperfine interaction. This internal 
field is a semiclassical approximation to the nuclear sys- 
tem, this approximation being valid in the limit of large 
number of nuclei, -/Vqd ^ 1. If the nuclei are noticeably 
polarized, this field almost coincides with the Overhauser 
field that represents the average nuclear polarization. It 
is important that the internal field persists even at zero 
polarization giving rise to doublet splitting of the order 

First we write the Hamiltonian in Eq. (^J) in the basis 
of the electron orbital states and present it as a sum of 
two terms 

-^hf = H HF + Vhf, (6) 
where the terms are defined as 

H° HF = AJ2\1)(1\S-J2\( R ^)\ 2 ^ (?) 

I k 

Vhf = ^lOa'l^-^ai^fcX^IO^. (8) 

l^l' k 

By definition H^ F does not couple different orbital levels. 
By combining Eqs. (|l|) and (0) one obtains the following 
Hamiltonian 



H = J2( £ i+WbB 1 -S)\1)(1\. 



(9) 
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We will regard the mixing term Vhf as a perturbation to 
H . The justification for this is that the typical fluctua- 
tions of the electron energy due to the hyperfine interac- 
tion are much smaller than the orbital energy separation, 
Ml. 

We now concentrate on Ho and formulate scmiclassical 
description of it. For this we consider the operator of the 
orbitally dependent effective magnetic field; 



where 



Bi = B + —Kt 

g^B 



(10) 



(ii) 



Our goal is to replace the operator Ki with a classical 
field. To prove the replacement is reasonable we calculate 
the average of the square for a given unpolarized nuclear 
state l/u) 



P\K?\») 



El 



QD 



-pf ] (1+1)// (12) 



where p\ = J drVq^Krll)] 271 is a dimensionless con- 
stant depending on the wave function of the orbital I. 
We cannot simply replace Ki by its eigenvalues: as in 
the case of the usual spin algebra, different components 
of Ki do not commute. In addition its square does not 
commute with individual components, [Kf, Kf] ^ 0. To 
estimate fluctuations of Ki we calculate the uncertainty 
relations between its components 



AK?AK? > 



N, 



3/2 
QD 



1 



N, 



1/2- 
QD 



(13) 



— 1/2 

Since Nq D < 1 we have proved that the quantum fluc- 
tuations in Ki are much smaller than its typical ampli- 
tude. The semiclassical picture introduced above is only 
valid for high temperatures, kT 3> EnN^^ where there 
are many states available to the nuclear system and the 

1 /2 

typical length is proportional to iVq D . For temperatures 
below E n NQ^ the nuclear system will predominantly be 
in the ground state and the classical picture breaks down. 
This is similar to the quantum mechanical description of 
a particle moving in a potential. At zero temperature it 
will be localized in some potential minimum and quan- 
tum mechanics will dominate. At sufficiently high tem- 
peratures the particle occupies higher energy states and 
its motion is well described by classical mechanics. Hav- 
ing established this we can replace the average over the 
density matrix of the nuclei by the average over a classi- 
cal field Ki, note that it has no 'hat', whose values are 
Gaussian distributed 



3/2 



expl 



3(K ; -K, (0) ) 2 

m?) 



) (14) 



Here is the average, or Overhauser, field assuming 

small degree of polarization. The effective magnetic field 
acting on the electron is Bi = Bin where n is the unit 
vector along the total field for a given Ki, see Fig. [l]. 
For this configuration the spin eigenfunctions are \n±) 
corresponding to eigenvalues 



n ■ S\n±) = ±^\ n ±) 
The effective Zeeman Hamiltonian is then 
Hz = gp-B Bi ■ S 



(15) 



(16) 



in this given internal field configuration. The spectrum 
of Hq thus consists of many doublets distinguished by 
the value of Bi . The magnitude of the effective field Bi 
determines the Zeeman splitting of the each doublet 



9PbBi 



2E B Kicost 



a/2 



(17) 



where Eg — gfisB the external field Zeeman energy. 

We conclude this section with two remarks concerning 
time and energy scales. First, since the dynamics of Ki is 
essentially the dynamics of the nuclear system it changes 
on a time .scale of the nuclear relaxation, which can be 
quite long.EEl For electron dynamics on shorter timescales 
Ki plays the role of a 'frozen disorder'. At longer time 
scales self-averaging over all configurations of Ki takes 
place. Second, the RMS value of Ki is approximately 
5 peV which corresponds to a 100 Gauss field (for a value 
N QD « 10 5 ). 



IV. TRANSITION RATE 

We concentrate on the transitions between the doublet 
components in Eq. (|l5| ) . Assuming that the higher energy 
doublet state is initially occupied, we will calculate the 
transition rate to the lower one. The transition must be 
accompanied by energy dissipation equal to Ao- piw, 
energy cannot be absorbed be the nuclear systemJiSEHl 
So we need an external mechanism of energy dissipation. 
The most effective one is known to be phonons. However 
the phonons alone cannot change the electron spin so we 
need a mechanism which mixes spin and orbital degrees 
of freedom, that is Vhf from Eq. (||) . Thus the transition 
amplitude is proportional to both Vhf and the electron 
phonon coupling -ff p h- 

Here we assume that the electron is in the lowest or- 
bital state |0) since the phonon mechanism will bring 
the electron to this state from any higher orbital on 
timescales much smaller than those related to transtions 
between the doublet components. Thus we consider an 
initial state of the entire system |i) = |0, n_; fj,; N) which 
is a product state of the electron, nuclear \u) and phonon 
\N) systems and the final state |f) = |0, n + ; /u'; AT'). The 
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transition amplitude between them, in second order per- 
turbation theory, reads 



r = E 

1^0 



(0, n+; h'\Vhf\1, n_; iV'|i? P h|0; JV) 



(s - El) + E B 
(OiJV^gphlii^^n+i^lVkFlO.n-iM) 

(£0 — Ei) — ^B 



(18) 



The summation is over virtual states involving higher 
orbitals and the denominators in Eq. (|is| ) contains the 
energy differences between different orbital states. The 
internal field depends on the orbital state, resulting in 
a rather complicated expression. Albeit the energy re- 
lated to the internal field is much smaller than the orbital 
separation so we can safely replace the Zeeman splitting 
with Eb- The reason for that is that only at high exter- 
nal fields where Aq ~ Eb will the effects of the Zeeman 



splitting be appreciable in the denominator. The internal 
field also appears in the phonon rate since it determines 
the electron energy difference between the initial and fi- 
nal state. Since -ff p h does not connect different nuclear 
states and conversely Vhf does not mix different phonon 
states the sums over intermediate phonon and nuclear 
states reduce to a single term. From this transition am- 
plitude the transition rate is obtained via Fermis golden 
rule 



(19) 



N> 



where E\ — Ef is the energy difference between the initial 
and final state of the combined systems. Substituting 
Eq. ( pi ) into Eq. ([ii]) we get the following relation for 
the spin-flip rate 



r sf = 



^ { ((so - iO + E B )((eo -e l>) + E B ) T £>' N ^ l > N)i ^ N ^ hl °' " ^ 

, (l,n.;^ F \0,n + )(l',n + \V HF \0,n.^) ^(0; N\H ph \l; N')(0; N'\H ph \l'-, N)S( Ei - El ) 



((e - ej) + E B )((e - e v ) - E B ) K 



N' 



(0 > n_; M |%F|i,tH>(0 I n + |% K |I' I n_;/*> N\H ph \Q; N')(l'; N'\H ph \0; N)6( Ei - E t ) 



((e -e t ) -E B )((s -e v ) + E B ) % 



N' 



(0,n_;/x|yHF|^n. + )(r,n + |yHF|0,n-;/x) 2tt ^ \ 

((£o _ £i) _ £ ; s ) ((eo _ ei , ) _ i , B ) T ^ NlHphl °> N){ °' N lHphl1 ' N)6(Ei ~ Ei) ) ' (20) 
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The spin-flip rate depends on the initial state of the nu- 
clear system |/x). This poses the problem of how to deal 
with the nuclear state since we already demoted all 
the spin operators to a collective classical variable. A 
conceptually simple solution lies in the fact that when 
Eq. ( pp| ) and (||) are considered together one sees that 
the rate is a sum over all pairs of nuclei in the system. 
Focusing on a given pair of nuclei k and k 1 , all the other 
nuclei are unchanged when the electron spin is 'scattered' 
on by this pair. By simply redefining the classical field 
such that it is composed of all nuclei except this given 
pair we can circumvent the problem. This procedure will 
not change our previous result regarding the properties of 



Ki and also defining |/x) = \mk)\mk>) makes it straight- 
forward to work with the nuclear states in Eq. (^0|) . 

Although the transition rate can be very slow, the typ- 
ical duration of a transition event is set by energy uncer- 
tainty, h/ Aq. This is much shorter than the typical time 
for nuclear spin relaxation so that the nuclear system is 
frozen in a given configuration of Kq. In this case taking 
an average over Kq is not well motivated since a given 
configuration of Kq is fixed. For now we will postpone 
the average over the classical field. After performing a 
thermal average over \fj,) in the transition rate we obtain 
the following 
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where Sei — Eq — £; . Here we have introduced the one-site 
nuclear spin correlation function G corr , which is evaluated 
in appendix |A| The denominator has been expanded to 
second order in the Zeeman splitting. The parameters 
aw, aw are related to the Uhf matrix elements 

aw = A 2 C n J d 3 R k *t(R k )\* (R k )\ 2 *v(R k ) (22) 

a w = A 2 C n J d 3 R k y*(R k )%(R k )V (R k ) 2 (23) 

and 7;;' , ^w are generalized phonon transition rates 

W - e^-- - 1 d(?kJ qi/ -Ao)(24) 

2^ ag(g)[e^-]o, ; [e^-]o,^ 
™' = T^- ifflfa^ - i S{Ku> qv -Ao)(25) 

qv 

Here we have only included the emission process since we 
assume the the spin is initially in the higher energy state. 

Until now we have considered a general quantum dot 
and the rate in Eq. ( pl| ) is valid for any quantum dot. To 
proceed further we will specify the confining potential 
to be parabolic in the lateral direction and in the trans- 
verse one a triangular well potential is chosen. The wave 



function is (r\l) = Xo{ z )*l J n,M{r,6) where n,M denote 
the orbital and angular momentum quantum numbers re- 
spectively of the Darwin- Fock solution and Xo( z ) is the 
wave function in the transverse direction. The square of 
the lateral confining length is I 2 — h 2 /m*hfl, where m* 
is the electron effective mass and fl — (Q 2 i + (oj c /2) 2 ) 1 / 2 is 
the effective confining frequency. The confining energy is 
Hflo and u) c = eB/m* is the cyclotron frequency. What 
remains is to calculate the a's and 7's in Eqs. (|22])-(p5|). 
The results of these calculations are presented appendix 

In principle it is possible to obtain the rate for all pa- 
rameter values but to make the discussion more transpar- 
ent we will consider regimes of the applies magnetic field 
(i) E B « E n N Q ^/ 2 and (m) E b > E n N Q ^ /2 . In regime 
(i) both Ao -C fic v i~ x and Ao -C kT (for experimentally 
relevant temperatures) and only the lowest order terms 
in Ao/hc^i^ 1 need to be considered. In GaAs ficil" 1 — 
3.3^ _1 nmxmeV and Ucti" 1 — 2.0^~ 1 nmxmeV for the 
longitudinal and transverse branches respectively. In the 
other regime we assume that the applied field dominates 
and that the internal field may be ignored, but no ad- 
ditional assumption are made in this case. This will 
in general require some numerical work. The resulting 
hyperfine-mediated spin-flip rates are 



^ - ^^ ^^ 2 ^^ {El + Kl + 2E B K QCaS 6) iorE B KE n N^ (26) 



HF 



G„„ E?(n(E B ) + 1) {eh u efE% /„,„,, (E B \* ,„ ,| . „ ^ „ „-i/2 



Note that the rates have different dependence on the 
emitted energy Ao. In regime (i) it is proportional to 
A 2 , due to the Bose statistics, and in regime (it) it is 
proportional to Ag for Ao <C hc^i^ 1 . In Eq. (|27| ) we 
introduce the functions Co and C2 which contain the de- 
tails of the higher orbitals and the anisotropy integrals. 
For low fields Ao <C hc^i -1 these functions are constant. 
The value for the spin flip rate in Eq. ( p6| ) , for some typi- 
cal value of K w 10~ 2 T, is very low r HF ~ 1CT 6 s. This 
results in a lifetime of days, which will be extremely dif- 
ficult to measure. We have plotted the general spin-flip 
rate in Eq. ( ^7| ) for different confining energies and tem- 
peratures in Figs. ^ and |j| for both in-plane and perpen- 
dicular applied magnetic field. In the case of the higher 
confining energy there is very small difference between 
the in-plane and perpendicular direction of the external 
magnetic field. For the lower confining energy there is a 
substantial difference between the two directions of mag- 
netic field. In this case the approximation Ao <C hc v l~ x 
is no longer good and the the Ag dependence of the rate 



is changed by the C-functions. The values of the rates 
are quite small, depending on the applied field, being in 
- Is- 1 for T = 4K at B » 0.5 T for a confining fre- 
quency of TiVLq = 2 meV. 



V. DISCUSSION 

Generally speaking spin-flip rates require external dis- 
sipation. This is why at small Zeeman splittings they 
will contain small factor reflecting the vanishing phonon 
density of states. For a spin-orbit rateE3, the Kramer's 
degeneracy results in this small factor being proportional 
to E g . Presence of nuclear spins violate the Kramer's 
theorem. So that the hyperfine rate discussed in the 
present paper is proportional to E B and will dominate 
at sufficiently low fields. 

Comparing the .-hjyperfine rate in Eq. ( p7| ) to spin- 
orbit related ratesi3 T^o we obtain that the hyperfine 
rate dominates for fields below B w 0.3 T (assuming 
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zq = 10 nm and Tl^Iq = 4meV). 

The role of the internal field produced by the nuclei is 
that the spin-flip rate does not vanish even for in the ab- 
sence of external magnetic field. The show that the min- 
imum rate is rather small, corresponding to a relaxation 
time of the order of days. We believe that the internal 
field will play an important role when the full dynamics 
of the electron spin in the presence of the nuclear system 
is considered. 

Our model should also be applicable to other polar 
semiconductors which have non-zero nuclear spin, e.g. 
InAs where the <?-factor is much larger. 
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= (S«_)*S p + _(6i a SP) T 6 k , k ,, (Al) 

where 5I a — I a — {I a )r and the electron spin matrix el- 
ements being S"_ = (n+\S a \n-). The kronekker delta 
reflects that there are no correlations between two dif- 
ferent nuclei and we have dropped the k subscript in 
(8I a 5I°)T since the nuclei is assumed to be identical. 
By defining the symmetric correlator 

g a P = ^(Si a 5i + 8i fi 8i a ) T (A2) 

we get the following 

Gcorr = {S%_Y Sl_{g aP + i/2e a ^(P) T ). (A3) 

In an isotropic system, (I)t — 0, the value of the corre- 
lation function is G corr = \ + 1). 



APPENDIX A: AVERAGE OVER NUCLEAR 
SPIN PAIR 

When the thermal average over a pair of free nuclei the 
following nuclear correlation function appear in Eq. (|21 



Gc 



(S?_)*S£_<Mj?M£) 



0' 



APPENDIX B: MATRIX ELEMENTS FOR A 
PARABOLIC QUANTUM DOT 

Using the Darwin- Fock solutions and the Fang- Howard 
variational solution for the triangular quantum well we 
obtain the following equation are Eqs. (p2|)-(^5|) 



. A 2 Cn(eh 14 l) 2 A 3 F{n + n' + \M\ + i)2-3(n+n'+|M|) 

VqdC 8?rpc 5 ^ n!n'!(n+ |M|)!(n' + \M\)\ 



( 



E 



id ( M 

cl \hc v 



2(n+n' + |M|-l) 



d{cosO)A v {6) 



( sin 0)2(n+n'+|M|) exp /_!(A^2\ 



(Bl) 



where C 1 = zq J cfe|x(.z)| 4 an( i c 5 = 



" t 5 is the 



effective sound velocity of the phonons. The anisotropy 
functions are 



MO) 
MO) 



9 cos 2 6 sin 4 9 



(B2) 
(B3) 



The equation for aw"fw is identical except for a dif- 
ferent Kronekker delta function 5m,— M'- The integral 
needs in general to be evaluated numerically but when 
A <C Tic u £~ l the exponential term and the denomina- 
tor become unity and the resulting integral is simple to 
calculate. 
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FIG. 1: The internal field coordinate system is set by the 
external magnetic field B, i.e. e z |] B. The combination of 
the external and the internal field Ki results in an effective 
field Bi. 
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FIG. 2: The hyperfine-mediated spin-flip rate for a quantum 
dot with zo = 10 nm and htlo = 5 meV, plotted a function 
of external magnetic field for two different temperatures T = 
0.1K and 4K. 
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FIG. 3: The hyperfine-mediated spin-flip rate for a quantum 
dot with zo = 10 nm and h£lo — 5 meV, plotted a function 
of external magnetic field for two different temperatures T — 
0.1 K and 4 K. 



